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GLOSSARY OF TERMS AND SYMBOLS 

Ai Incremental area of one and two dimensional 

velocity distribution grids 

B Slot jet width 

S Distance between jet centerlines (units in 

feet) 

D Jet diameter (ft.) 

Q Volume flux (cu.ft./sec) 

Qo Initial volume flux (cu.ft./sec) 

J Momentum flux (ft. ^ /sec 2 ) 

A 2 

Jo Inital momentum flux (ft. /sec ) 

X Axial distance downstream of jet source. See 

Fig. 1 

Y Direction parallel to axis perpendicular to 

and through centerlines of the row of jets. 

See fig. 1 

Z Direction perpendicular to X-Y plane of jet 

array. See fig. 1 

Jet array Row of jets 

Radial Perpendicular to jet axis 

U Velocity component along jet axis or X axis 

(ft ./sec) 

V Velocity component along Y axis 

W Velocity component along Z axis 

Uo Maximum initial jet velocity as measured at 
jet source (ft. /sec) 

Urn Maximum axial jet velocity as measured at a 
given downstream cross-section (ft. /sec) 

U(r=l/2Ro) Axial velocity measured at midpoint bet- 
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ween jet centerline and nozzle wall at 
jet source. 

Section Y-Z plane in which velocity distributions 

fall 

Ui Average of 200 instantaneous velocity 

readings that is assigned to a one or two 
dimensional coordinate in the velocity distri- 
bution grid (ft. /sec) 

Uiave Average velocity found by [U ( i ) + U(i+l)]/2 
in the one dimensional grid and by [U(i,j) + 
U(i+l,j) +U(i,j+1) + U (i+1, j+1) ]/4 for the 
two dimensional velocity distribution grid, 
(ft ./sec) 

Re Reynold's number, based on maximum nozzle 

velocity at the jet source and on the jet 
diameter 

r Correlation coefficient 
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Chapter I 
Introduction 


The subject of this study is the free jet which 
results from the injection of a fluid through a row of 
equally spaced nozzles (at turbulent Reynold's numbers) 
into an unbounded resevoir in which a second fluid is 
initially at rest. The focus of attention is on one 
particular property of the free jet - the mass flow 
rate across a section at right angles to the jet axis: 
m = J~ p U dA eqn. (1) 

As fluid is injected into the quiescent environ- 
ment, the second fluid of this resevoir is drawn 
radially inwards towards the jet axis (in the case of 
the circular jet) or is drawn in along the Z axis (in 
the case of the slot jet) until it becomes part of the 
jet. This process is known as entrainment. 

The entrainment rate for a fluid issuing into a 
quiescent environment can be varied between the limi- 
ting entrainment rates of the single circular and slot 
jet by changing the jet spacing to jet diameter ratio 
(S/D). Choosing smaller S/D's results in entrainment 
rates closer to that of the slot jet, while larger 
S/D's result in entrainment rates which approach that 
of the circular jet. 
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In a study of dilution jet mixing in compact 
reverse flow combustors for gas turbines, Lipshitz and 
Greber noted the need for and absence of information on 
the entrainment rate for a row of jets (ref. 6). 

It is the intent of this study to provide funda- 
mental data and its results on the effect that jet 
spacing in a row of jets has upon the entrainment rate. 
Chimney spacing required to givea specified dilution 
rate of exhaust gases into the environment, fuel to air 
mixing ratios in injection and carburation systems, the 
determination of jet interaction behavior in cases such 
as V/STOL aircraft and jet engine combustion chambers 
represent only a few of the many applications that 
require basic information on jet array entrainment 
rates . 

This study undertakes as its primary objective the 
establishment of entrainment rates for array spacing 
ratios of S/D = 2.5, 5, 10, and 20. Velocity distribu- 
tions are determined at X/D = 16.4, 32.8, 65.6, 98.4, 
131.2, and 164 at Reynold's numbers ranging from 5,110 
to 17,000. Centerline velocity decay rates are also 
presented in the results. 



Chapter II 
PREVIOUS STUDIES 


Numerous measurements of the axial velocity pro- 
files have been made in studies of isothermal air jets. 
Inserting U(y,z) into equation (1), the mass flow can 
be determined by numerical quadrature. This is, how- 
ever, a time consuming and difficult procedure that 
many authors have circumvented by assuming a profile 
that permits analytical integration. Then by measuring 
the jet velocity at the axis and the distance from the 
axis at which the velocity is half of the center line 
velocity, the assumed velocity profile can then be 
scaled to that section. Obviously, the mass flux then 
becomes dependent on the 'theoretical' curve chosen by 
the author. 

The difficulties in determining detailed velocity 
distributions are largely a result of the limitations 
of the instrument used to make the velocity measure- 
ments - the pitot tube (ref. 9). At large distances 
from the jet axis, where the velocity of jet becomes 
very small and increasingly intermittant, the pressure 
differences are difficult to read because of their 
small and constantly changing magnitude. 

The limitations involved in determining detailed 
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velocity distributions and in curve fitting the veloci- 
ty profile into an analytical form that can than be 
integrated prompted Ricou and Spalding to formulate an 
experimental procedure that bypasses the integration of 
a velocity profile altogether. The two researchers 
devised their method on the basis of the following 
considerations (ref. 9): 

"In the ideal free jet which is under consid- 
eration, the surrounding resevoir is large and, 
except in the vicinity of the jet axis, at uni- 
form pressure; the entrained fluid flows radially 
inward towards the jet axis. If, however, a tur- 
bulent jet is partially enclosed in circular walls, 
the radial inflow is impeded; the fluid entrained 
by the enclosed part of the jet has to flow axial- 
ly in the annular space between the cylindrical 
wall and the conical 'boundary' of the jet; the 
corresponding axial pressure gradients can be 
be measured. 

Now let us suppose that the cylindrical wall is 
made porous, and that a controlled and measured 
amount of fluid can be caused to flow through it 
in a radially inward direction. Then, if this flow 
is equal in quantity to that which would have 
passed through the area occupied by the wall if 
the wall had been absent, the axial pressure 
gradients reffered to in the last paragraph would 
disappear . 

This recognition underlies the experimental 
method which has been used: the flow rate through 
the porous wall is varied until no axial pressure 
gradients can be detected; this flow rate is then 
measured and presumed to be equal to that which 
would be entrained if the jet were unenclosed." 


A detailed discription of their apparatus and procedure 
may be found in reference 9. 

If we take the case of the isothermal free jet - where 
the temperature of the resevoir air is equal to that of 
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fluid issuing into the resevoir - the entrainment rate 
reduces to the rate of change in the volume flux. The 
volume flux itself can be expressed as a function of 
the initial volume flux, Qo, the distance from the jet 
source, X, the jet diameter or width, D or B, and the 
jet Reynold's number. Re, and the initial velocity 
profile. In nondimensional form, the expression becomes: 
Q/Qo = f (X/D or X/B, Re, initial profile) 

For a given nozzle geometry, a measure of the 
initial velocity profile, called the coefficient of jet 
structure, has been defined as a = 0.066(Um/Uo) where 
Um/Uo is a centerline velocity divided by the maximum 
initial velocity in a jet whose initial velocity pro- 
file is not uniform (fig. 2). The coefficient of jet 
structure increases from a base value of 0.066 (when 
the initial velocity profile is uniform) with increas- 
ing nonuniformity in the initial velocity profile. The 
product of a/0.066 and Um/Uo (where Urn is a centerline 
velocity for a jet whose initial velocity profile is 
not uniform) now agrees with the theoretical curve (due 
to Tollmien) for the variation of velocity on the jet 
centerline. The dependence of the coefficient of jet 
structure on the initial velocity profile is shown in 
figure 3. The volume flux can then be expressed as: 

Q/Qo = f (X/D or X/B, initial velocity profile) 
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Figure 25 shows a centerline velocity, decay rate 
curve for the circular jet that has been adjusted to 
the initial velocity profile of the S/D = 20 curve just 
above. The cor respondance was so exact that the circu- 
lar jet curve had to be shifted slightly downward in 
order to distinguish the two curves. 

Analysis of the experiments of Trupel, Zimm, the 
Gottingham Aerodynamics Institute, Turkus, and Syrkin 
(ref. 1) has shown that for Reynold's numbers ranging 
from 20,000 (or where the initial velocity distribution 
becomes uniform) to 4,000,000 the plot of Um/Uo versus 
X/D falls on the same curve, and the entrainment rate 
can be reduced to: 

Q/Qc = f (X/D or X/B) 

As the Reynold's numbers decrease and the nonuniformity 
in the initial velocity distribution increases, the 
axial velocities downstream of the jet source decay 
more rapidly. The result is that the volume flux rates 
of turbulent jets with an Re of less than 20,000 are 
smaller than entrainment rates of turbulent jets with 
uniform initial profiles. It should be noted that the 
formulas given for the slot and single circular jet 
(table 1) are generally valid for Reynold's numbers 
in excess of 20,000. 

The entrainment rates for the row of jets requires the 
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addition of the nondimensional parameter that takes 
into account the jet spacing. The functional relation 
for the row of jets is therefore represented by: 

Q/Qo = f (X/D, S/D, initial velocity profile) 

Mary Nash, a graduate student at Case Institute of 
Technology, presented results for two spacing configu- 
rations (S/D = 5 and 10) for the row of jets. Her 
results are: 

Q/Qo = (X/D)' 665 
for S/D = 5, at Re = 5600, and 

Q/Qo = (X/D)' 692 

for S/D = 10, at Re = 7,000. These entrainment rates 
are limited to an X/D range of 70 to 220. For X/D < 70, 
Nash found her experimental values for volume fluxes 
(for both S/D = 5 and 10) to be higher than the rate 
for the single circular jet. As a result, she deemed 
her data valid only for X/D > 70. 


Chapter III 
EXPERIMENTAL METHOD 

The volume flux is determined at the jet source, 
Qo, and at cross sections 5, 10, 20, 30, 40, and 50 cm. 
from the jet source. The volume flux, £ &Y 2 Uiave A Z, 
is found by integrating the velocity distributions that 
have been measured across the entire width of the jet 
and from the midpoint between the central jet being 
surveyed and its neighbor underneath to the midpoint 
between the center jet and the neighboring jet above, 
as shown in figure 4. In all cases, the hot wire probe 
is oriented along the Y axis. 


2.1 Apparatus 

41 jet apertures, .3 inches on center, are drilled 
and reamed to .12 inches in diameter in a 2.25 x 13 x 1 
inch block of polyvinylchloride. Thisspacing between 
adjacent jets gives an original S/D ratio of 2.5, with 
larger S/D ratios obtained by simply plugging ever 
other jet. The row of jets is mounted on an adjustable 
support, and is connected to a plenum chamber that 
evenly divides the flow from the compressor (supplied 
through 1 inch O.D. tubing) to each one of the aper- 
tures . 

The row of jets is bounded by a three sided enclo- 
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sure (figure 5) that constrains the entrainment air to 
be drawn in along the Z axis, and that prevents the 
flow from spreading along the Y axis. The setup not 
only approximates flow in a section of an infinite row 
of jets, but prohibits the axial entrainment of air 
from behind the jet source. The enclosure is construct- 
ed out of quarter inch plexiglas, glued together with 
thermal adhesive. 

A three foot diameter axial compressor supplies 
air to the plenum chamber through a Jamesbury valve 
used to adjust flow as read on a 28 cfm. full scale, 
Fischer-Porter rotameter. 

A three foot hot wire probe support is mounted on 
a three directional tracking mechanism that is powered 
along the Z and Y axis by Superior Electric MO62-FD03 
stepping motors. The stepping motors are connected to 
screw and follower mechanisms that drive the probe at 1 
inch per 16 shaft revolutions along the Z axis and at 1 
inch per 13 shaft revolutions along the Y axis. Opera- 
ting at 1.8 degree steps per pulse, the steppers are 
capable of moving the probe incremental distances which 
are on the order of the thickness of the hot wire probe 
diameter - .002 inches. Even at such small distances, 
the error due to hysterisis in the follower mechanism, 
a double follower is employed in the Z direction 
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mechanism/ while the weight of the probe mount assures 
no hysterisis along the Y axis (figure 6). Smooth 
tracking/ resulting in minimal torque required from 
the stepping motors, is obtained by mounting the Z and 
Y axis traversing components on linear tracking bearings. 

A Tecmar Labmaster A/D board supplies the square 
wave pulses that are amplified by a single driver card 
which powers the stepper motors. A pulse rate of just 
under 400 per second (the maximum allowable for accu- 
rate operation of the stepping motors) is generated by 
a timer source on the A to D board, and is in turn 
gated by a second timer that determines the duration of 
the square wave form. In this way, the distance that 
the probe is moved can be controlled by varying the 
duration set by the second or gating timer. 

The probe support is fixed to the traversing 
mechanism so that it moves the hot wire probe into the 
jet wake along the Z axis. This probe support orienta- 
tion places the traversing mechanism the farthest pos- 
sible distance from the jet, while minimising the ob- 
struction that the probe support itself introduces into 
the flow field of the jet. 

2.2 Instrumentation 

A TSI model 1054B constant temperature linearized 
hot wire anemometer is paired with a model 1210-20 hot 
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wire probe that produces a voltage signal which is 
proportional to air speed. The analog voltage signal 
generated by the anemometer is then sampled by the 
Labmaster A/D board at a 50 hertz frequency, and pro- 
cessed by an IBM XT personal computer (256k with one 
hard and one floppy disk drive) . 

All plots are made on a Hewlett-Packard 7470A 
plotter . 
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2.3 Experimental Procedure 

The initial volume and momentum fluxes are experi- 
mentally determined at the exit of the jet nozzle 
through the establishment of a highly detailed velocity 
distribution across the jet nozzle. The hot wire probe 
sensor element is positioned at centerline of the jet, 
just inside the nozzle. The program is initialized, and 
the probe is positioned at .003 inch increments across 
the mouth along the Z axis (from the jet centerline to 
the nozzle wall) of the jet for velocity readings. The 
values for Co and Jo are calculated by the following 
formulas : 

Qo = X (Uiave) (Ai) 

Jo = £ (Uiave) 2 (Ai) 

Ai = 2 [R Ci ) 1 ~ R(i-l) 2 J 

The volume flux found by integrating the initial velo- 
city profile is compared to the volume flux as indi- 
cated by the rotameter. 

Downstream cross sections at 5, 10, 20, 30, 40, 
and 50 cm. are surveyed by locating the jet ' boundary' 
using smoke to visualize the flow, positioning the hot 
wire probe just outside the jet boundary, and initiali- 
zing a program that positions the probe at .0625 inch 
intervals across the jet wake in the Z direction until 
the probe is 1/4 inch outside the jet 'boundary' on the 
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otherside. The position of the probe is incremented up 
in the Y direction by .075 inches and the Z direction 
operation is performed again until the the probe has 
been moved up from the midpoint between the central jet 
and it neighbor below to the midpoint between the 
central jet and its neighbor above. The probe therefore 
sweeps out the area as shown in figure 4: where the 
area equals the distance from boundary to boundary 
times S. This procedure is repeated for all cross 
sections, although in S/D ratios of 10 and 20 the 
velocity grid is relaxed along the Y axis to .15 inches 
at far downstream cross sections. If, at any point in 
this operation, there is a question about the exact lo- 
cation of the jet 'boundary', the program is momentari- 
ly stopped and the flow field of the jet is again 
visualized using smoke. 

Having established values for the volume flux at 
the jet source and at the six downstream locations, 
functions of the form 

Q/Qo = A (X/D) c + B 

are determined by choosing the exponent, C, and then 
performing a least squares fit which determines the 
constants A and B. 



Chapter IV 
RESULTS 

The entrainment rates found for the four spacing 
ratios tested are presented in Table 2, along with the 
previously determined entrainment rates for the slot 
jet and the circular jet. Figure 7 is a graphical 
representation of the rates given in Table 2. Figures 8 
- 11 show the least squares fit of the entrainment rate 
curves to the experimentally determined volume fluxes 
for each of the spacing ratios. As expected, the en- 
trainment rates increase monot onicai ly with increasing 
S/D, with small S/D jet rows approaching the entrain- 
ment behavior of the slot jet and large S/D jet rows 
approching the entrainment behavior of the circular 
jet . 

The entrainment rates for the four spacing ratios 
are valid for X/D greater than 10, because self-simi- 
larity in the velocity distributions generally begins 
at 10 nozzle diameters from the jet source (ref. 1). 
For S/D = 5, 10, and 20, the entrainment rates between 
X/D = 10 and lower value delineating the valid range 
for the use of the given rate formula may be calculated 
by using the entrainment rate formula for the circular 
jet . 

For each of the four spacing to nozzle diameter 
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ratios tested, the upstream entrainment rate approaches 
that of the individual circular jet, while the down- 
stream entrainment rate approximates that of the slot 
jet . 

The least squares fit formulas for the velocity 
decay rates are presented in table 3, and are graph- 
ically represented in figure 12. The velocity decay 
rate formulas for each of the spacing ratios are graph- 
ed against the measured velocities in figures 13 - 16. 
The decay rates also vary monotonical ly with 
increasing S/D between the limiting cases of the slot 
and circular jet, with the decay rate of the smaller 
S/D row of jets approaching that of the slot jet and 
the decay rate of the larger S/D row of jets approach- 
ing the circular jet. 

The dimensionless plots of Q/Qo versus Z/Z(l/2Um) 
in figures 17 - 20 exhibit the se 1 f -simi 1 ar ity of the 
velocity distributions recorded for each of four spa- 
cing ratios. 

Summaries of directly measured and calculated 
quantities - listing values for S, Qo, Jo, Uo, Re, Urn, 
J, and Q - are given in tables 4-7. 

Figure 21 presents a plot of Q/Qo (experimentally 
determined) versus S/D for X/D = 65.6, 98.4, 131.2, and 
164. A least squares curve is fit to the dimensionless 
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volume fluxes for each X/D in order to predict dimen- 
sionless volume fluxes for S/D's ranging from 2.5 to 
20. The X/D = 164 curve is proportional to (S/D)* 65 , 
the X/D = 131.2 curve is proportional to (S/D)* 60 , the 
X/D = 98.4 curve is proportional to (S/D)* 55 , while the 
X/D = 65.6 curve is proportional to the square root of 
S/D. The curves for X/D = 164 and 131.2 are proportion- 
al to S/D in the S/D range between 2.5 and 10. 

I 

I 



CHAPTER V 


DISCUSSION 

The early stages of the present study were spent 
in an effort to assure the proper alignment of the jets 
first with respect to each other and second with res- 
pect to the traversing mechanism that moves the hot 
wire probe along each of the three axes. Misalignment 
in the row of jets itself can cause significant varia- 
tions in centerline velocity and volume flux values for 
a section. These values will be larger than those for 
the perfectly aligned case if the nozzles converge in 
the X - Y plane, while the centerline velocities and 
volume fluxes for the divergent nozzles will be smaller. 

Problems generated by jet misalignment were so 
severe that an original jet row assembly consisting of 
drilled and reamed threaded rod that screwed into tap- 
ped aluminum plate was discarded in favor of a solid 
block of aluminum in which jet apertures were drilled 
and reamed on a vertical milling machine. A number of 
jet row assemblies were drilled and tested, each one 
showing better symmetry in recorded velocities about 
the jet centerline along the Z axis and more uniform 
periodicity in velocities measured through the jet 
centerlines along the Y axis. Misalignment was further 
reduced by using a block of polyvinylchloride. The 
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softer plastic reduced drill over heating due to fric- 
tion and required less load on the drill thereby redu- 
cing drill deflection during the machining operation. 

After aligning the final row of jets assembly with 
traversing mechanism tests were performed to determine 
symmetry in velocity distributions about the jet cen- 
terline and to determine periodicity of velocity mea- 
surements taken through the jet centerlines along the Y 
axis. Figure 22 plots velocity measurements taken along 
the Y axis (through the jet centerlines) against the 
jet centerlines - as represented by the lines along the 
ordinate. The S/D = 5 case (shown in the lower half of 
the graph) shows the velocity distribution - at X/D = 

14 - for 13 central jets out of a 20 jet row. The S/D 
= 10 case represents the velocity distribution - at X/D 
= 65.6 - for the 7 central jets in a 12 jet row. Tests 
conducted at X/D = 164 for S/D = 2.5, 5, and 10 resul- 
ted in velocity distributions (along the Y axis) that 
are uniform throughout a central 11 inch portion of the 
13 inch distance from the top to bottom wall of the row 
of jet enclosure. 

There is a fundamental difference in the velocity 
profiles of the slot and circular jet at and beyond 
Y/Y(l/2Um) = ± 1.85 suggested in the comparison of 
figures 23 and 24. In both cases, the theoretical curve 
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due to Tollmien predicts a gradual decrease in axial 
velocity from Y/Y(l/2Um) = 2.0 until it reaches a value 
of zero at Y/Y(l/2Um) = 2.5. For the circular jet, ex- 
perimental data is in excellent agreement with 
Tollmien. Experimental data on the slot jet, however, 
indicates that between Y/Y(l/2Um) = 1.6 and 1.85 axial 
velocities are presented only to .2 to .l(Um), after 
which no data appears to exist. Possible reasons for 
this will be discussed below. Figure 25 shows qualita- 
tively the streamline patterns for a circular, free 
turbulent jet. At the jet 'boundary', marked by the 
broken line, the streamlines vary from nearly radial to 
nearly axial in direction. Farther downstream, this 
streamline direction change is more gradual. There 
seems to be little documentation in the literature on 
streamline behavior for the slot jet. Figure 26 illus- 
trates the results of flow visualization (using smoke) 
for row of jet spacing to diameter ratios of up to and 
including S/D = 10, and strongly suggests - as does 
figure 23 - that there exists a point where the axial 
velocity distribution goes to zero. 

The flow visualization used in this study indi- 
cates that there is not a monotonic decrease in axial 
velocities to zero (as in the case of the circular 
jet), but rather a region of oscillatory flow, as well 
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as back flow at a fairly well identified jet 'boun- 
dary'. Furthermore/ it has been found that there is a 
correlation between the 'boundary' and the hot wire 
anemometer voltage response at that location in the 
flow field. Looking again at figure 26, the mean 
streamlines approaching and crossing the jet 'boundary' 
(delimited by the broken line) are represented by a 
smooth lines v*ith both a radial and a somewhat smaller 
negative axial component. As this entrained air ap- 
proaches the jet 'boundary', the back flow intensifies 
until it becomes parallel to the 'boundary'. The back 
flow then ceases - as shown in figure 26 - oscillates 
back and forth along the X axis with no net axial 
velocity and with only a very small radial velocity. At 
this point, Ui recordings drop from axial velocities of 
approximately 2 fps. to readings that fluctuate between 
1.5 to .8 fps.. This represents the zone just outside 
the jet 'boundary' where no net axial velocity exists. 
Because the anemometer voltage signal records only the 
magnitude of the velocity and not the direction, an ap- 
parent net velocity is indicated when in fact none 
exists. As the probe moves farther out from the jet 
axis, Ui readings rise again to 1.2 to 1.8 fps.. These 
slightly higher readings show no fluctuation and repre- 
sent the zone in which the back flow' of the entrainment 
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air suddenly intensifies. Beyond this zone, velocity 
readings fall back to under 1.0 fps. from which they 
slowly decrease as the probe continues to move out from 
the jet axis. 

The velocity distributions shown in figures IB - 
52B represent total velocities. Within what has been 
termed the jet 'boundary 1 , flow visualization has shown 
that the velocities are axial. For the X/D locations of 
65.6 only velocity points within Y/Y(l/2Um) = ± 1.85 
are summed over their respective area to yield a volume 
flux. Farther downstream, the Y/Y(l/2Um) bounds shrink, 
until at X/D = 164 only the velocity points within 
Y/Y(l/2Um) = ± 1.6 are included in the axial velocity 
distribution. Since the downstream flow behavior of the 
row of jets approaches that of the slot jet, it is 
encouraging to note the very close cor respondance bet- 
ween the velocity distribution cutoffs - as identified 
by the flow visualization used in this study - and the 
experimental data on the slot jet from previous re- 
search. In the same way, the cutoff points for the X/D 
= 16.4 and 32.8 cross sections are approximately at 
Y/Y(l/2Um) = 2.6 and at Y/Y(l/2Um) = 2.2 respective- 

ly. In the case of the more closely packed row of jets 
- where streamline behavior approaches that of the slot 
jet - it is easier to identify the characteristic sig- 
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nature recorded by the anemometer because the back flow 
and oscillating flow at the 'boundary' are more pro- 
nounced than in the case of the larger S/D row of jets, 
whose streamline behavior approaches that of the circu- 
lar jet. 

Figures 27 - 35 display the flow behavior of en- 
trained air as it is drawn in towards the jet axes. In 
each of the figures, the row of jets is on the far 
right hand side of the photographs. The line of sight 
is almost at a right angle to the jet axis, and about 6 
inches below the level at which the smoke is being 
injected. The locations of the smoke tests in the jet 
flow field are indicated in figure 26: for the X/D = 30 
locations, smoke is injected into the jet as marked by 
"1", and for X/D = 120 as marked by "2". Only for the 
circular jet (figure 25) is the back flow - that chara- 
cterizes the entrainment flow for the row of jets - 
absent. Even as far upstream as X/D = 32.8, the S/D = 
20 row of jets exhibits a very slight back flow, while 
farther downstream (approximately X/D = 120) where sig- 
nificant mixing between the jets has begun the entrain- 
ment air is slowly oscillating along the X axis with a 
small net back flow. As the S/D ratio decreases, the 
strength of the back flew - both near upstream and 
downstream - increases. Compilation of data on the 
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initial momentum flux and its values for downstream 
cross sections of a row of jets suggests instead that, 
because very near upstream values and far downstream 
values of J can vary by as much as 5 - 10% from Jo, 
adjusting the volume flux to correspond with distri- 
bution that yield agreement with Jo would seriously 
inflate the near upstream and far downstream values of 
Q. Furthermore, the existence of backward entrainment 
flow along the boundary of the jet would indicate that 
the assumption upon which the conservation of axial 
momentum is founded - i.e. that the pressure everywhere 
in the flow field is uniform and ambient - is not 
entirely true. This study strongly suggests that pres- 
sure driven flows are in evidence at the boundary of 
the row of jets, and that consequently the conservation 
of axial momentum cannot be relied upon to furnish a 
cutoff for the axial velocity profile. 

At S/D = 2.5 and 5 tests, the initial volume flux 
as calculated by integrating the initial velocity pro- 
file corresponded to the volume flux as read by the 
rotameter to within 2%. This excellent agreement broke 
down in the S/D = 10 and 20 tests when higher plenum 
chamber pressures apparently resulted in leaks that 
caused rotometer readings to indicate volume fluxes up 
to 10% higher than those calculated by integrating the 
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initial velocity profile. No other explanation is read- 
ily available, sinse calibration of the rotameter indi- 
cated perfect linearity in its 0 to 28 cfm. readings. 
Because of the possibility of small leaks and the 6% 
inaccuracy associated with reading the rotameter, the 
initial volume flux as calculated by integrating the 
initial velocity was relied upon as primary informa- 
tion . 

The results (given in table 2) indicate that the 
entrainment for the row of jets behaves like the circu- 
lar jet in the near upstream, while behaving like the 
slot jet downstream. The 'transition' point, marked by 
the conjunction of the circular jet entrainment curve 
and the curves that are proportional to the square root 
of X/D, corresponds to an X/D at which the jets have 
mixed sufficiently so that all velocities along the 
Y axis, through the jet centerlines at that section, 
are axial. Using data on the jet width growth from the 
circular jet tests to predict the X/D at which this jet 
mixing would take place for each of the spacing ratios, 
excellent correspondance was found between the experi- 
mental and predicted value of X/D at which the jet 
'boundaries' cross. The predicted value is within 10% 
of the experimental value. 


CHAPTER VI 


RECOMMENDATIONS 

The following general recommendations concerning 
experiments on arrays of jets are based on the exper- 
ience resulting from the experiments of this study. 

The use of a three sensor probe anemometer would 
facilitate velocity measurement sinse it records the 
instantaneous direction as well as the speed of a flow. 

Flow visualization techniques should be employed 
from the outset of the study so that velocity measure- 
ments may be interpreted in light of a previous under- 
standing of general flow characteristics. 

Extreme care must be taken in aligning the jets 
with respect to each other and with respect to the 
traversing mechanism. 

Finally, it is useful to work with initial velo- 
city profiles sinse entrainment rates found will gene- 
raly be valid over a wide range of Reynold's numbers in 
which the initial velocity profile is uniform. 
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TABLE 1> 


Turbulent 

otudioo 


jot rotoo from proviouo 
<rm f . 2. 7. and 9) 


Jot Typo 

Entrainmont 

rato 

Contorl ino 
Volocity 

Slot Jot 

Q/Qo - . 62 \Jx7b 

Um/Uo - 2. 28'J§7X 

Circular Jot 

0/Qo - . 32 <X/D> 

Um/Uo - 6. 2CD/X) 

Row of Joto 
S/D - 5 

o/oo - cx/oy 665 

Um/Uo - . 24 CD/X1* 17 

Row of Joto 
S/D - 10 

gQ9 

Q/Oo - <X/0> 

Um/Uo - . 74 <0/X)- 50 
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TABLE 2s Erstrainunt Rotu 
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TABLE 3« Cantor lino Velocity Rates 


S/D 

Re 

Velocity Decay Rate 

Range 

(X/D) 

1 

2.5 

5170 

Um/Uo - 1.53C0/X)* S5 +. 0030 

16. 4-164 

. 999 

5 

8800 

Um/Uo - 3. 65<D/Xr 90 0169 

16. 4-66 

.998 

5 

8800 

Um/Uo - 1. 36<D/X5* 60 -. 0068 

32. 8-164 

.999 

10 

9730 

Um/Uo - 5. 21 <D/X)' 95 0079 

16. 4-66 

. 998 

10 

9730 

Um/Uo - 2. 39<D/XX 80 +. 0060 

65. 6-164 

. 998 

20 

12070 

Um/Uo - 4. 77 <D/X>* 95 0030 

16. 4-164 

. 999 
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TABLE 4i Summary af Data for S/D - 2. 5 


S/D -2.5 


s 

Qo 

Jo 

Uo 

Ra 

. 025 

.0057 

. 425 

84. 5 

5170 

X/D 

Urn 

J 

Q 





16. 4 

28. 07 

. 429 

.023 

32.8 

18. 79 

. 448 

.035 

65.6 

12. 58 

. 430 

.050 

90. 4 

10. 14 

.414 

.058 

131.2 

9. 32 

. 434 

. 067 






164 


7. 92 


417 


074 
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TABLE 5i Summary of Data for S/D ■ 5 


] ! 

I S/D - 5 


i i 

Qo 

Jo 

Uo 

Rq j 

• 




.050 1 

i 1 

. DID 



1. 07 

142. 5 

1 

8800 ! 

i 


j 

X/D 

Um 

j 

! 

o ! 

i 





16. 4 

. 

44. 81 

.995 

. 054 

32. 8 

» 

23. 86 

1.06 

.073 | 

t 

' 

65. 6 

_ _ _ 

15. 20 

1. 16 

1 

. 113 ! 

i 

"■ 

98. 4 

12. 24 

1. 07 

1 

. 133 

i 

131.2 ! 10.07 

1 

1.04 | .147 

i J 

164 j 8. 98 

i — . L 

1. 00 

1 

. 163 j 

J 
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TABLE 61 S ummar y of Data for S/D ” 10 


S/D - 10 


Qo 


Jo 


llo 


Rq 


. 100 


0110 


1.30 


158.0 


9730 


1 1 1 ■» ' ■ - - 

X/D 

Urn 

J 

Q 

16.4 

57.00 

1.30 

.057 

32.8 

28.40 

1.29 

. Ill 

65.6 

14.33 

1.30 

. 159 

98. 4 

10. 36 

1. 39 

.203 

131. 2 

8. 40 

1. 30 

.22 6 

164 

7. 52 

1. 30 

.253 
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TABLE 7s Summary of Data for S/D m 20 

S/D - 2D 

S Qo Jo Uo Rq 

.200 .0132 2.16 1 96.5 12070 j 

J i L____ 

X/D Urn J • Q 

16. 4 64.89 1.90 .068 

j 

32.8 33.44 2.20 -135 i 

i 

1 4 j 

65.6 17.43 j 2.19 .255 i 
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FIGURE 2i Variation of velocity on the axle 
of a circular Jet (ref. 1) 
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FIGURE 3i The dependence of the coefficient of 
structure of o Jet, a, on the initial 
velocity profile (ref. 1) 
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ORIGINAL' PAGE 

BLACK AND WHITE PHOTOGRAPH 


FIGURE 5: Row of jets, enclosure, and smoke nozzle 
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ORIGINAL PAGE! 

BLACK AND WHITE PHOTOGRAPH 



FIGURE 6* Row of Jots, enclosure, and 
traversing mechanism 
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FIGURE 12« Velocity dacay rata* For diffarant 
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FIGURE 23* Velocity distribution in a two-dimeneion- 
al turbulent Jet. Measurements due to Foerthmann 
Curve 1 due to theory of To 11 mien (ref. 10 . p. 747) 
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FIGURE 24* Velocity distribution in a circular, tur- 
bulent Jet. Measur e ments due to Relehardt. Curve 
1 due to theory of Tollmien (ref. 10. p. 749) 
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FIGURE 
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ORIGINAL PAGE 

BLACK AND WHITE PHOTOGRAPH 



FIGURE 28: Smoka tost, for S/D ■ 20 

at X/D » 30 
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ORIGINAL' PAGE 


BLACK AND WhT 
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FIGURE 29* Smoko tost for S/D ■ 20 

at X/D - 120 
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original page 

Black and white photograph 



FIGURE 30* Smoko tost for S/O ■ 10 

at X/D - 30 
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FIGURE 31 i Smoko toot for S/0 ■ 10 

at X/D - 120 
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BLACK AND WHITE PHOTOGRAPH 



FIGURE 32* Smoko tost for S/D ■ 5 

at X/D - 30 
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BLACK and WHITE photograph 



FIGURE 33* Smoke test for S/D - 5 

at X/D - 120 
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FIGURE 34, Smoko tost for S/D -2.5 

at X/D - 30 


65 


ORIGINAL PAGE 
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FIGURE 35* Smoke test for S/D ■* 2. 5 

at X/D - 120 
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APPENDIX A 

Error Due to Experimental Apparatus 

The early stages of the present study were spent 
in an attempt to minimize the effects of misalignment 
of the jet nozz 1 es first with respect to each other and 
second with respect to the traversing mechanism which 
moves the “Tiot wire probe along each of the three axes. 
Misalignment in the row of jets itself results in 
momentum losses that - while difficult to quantify - 
are undoubtedly very significant. Because of the pro- 
blems generated by jet misalignment, a drilled, 
threaded rod and tapped aluminum plate assembly was 
discarded in favor of a solid block of polyvinylchlo- 
ride in which jet apertures were drilled and reamed on 
a vertical milling machine. 

The length of the row of jets introduces another 
problem in so far as the study is ideally one of an 
infinite row of jets. Limitations in compressor flow 
rates in turn limit the length of the jet row, sinse a 
minimum turbulent Reynold's number must be maintained. 
Consequently, barriers were introduced at either end of 
the row of jets to prevent spreading along the Y axis. 
The unfortunate side effect of the enclosure is that it 
introduces a shearing force at the walls. In an attempt 
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to detect any influence that the walls might have on 
the central portion of the flow, the hot wire probe was 
stationed at a location far enough downstream that the 
flow was essentially two dimensional. The probe was 
then traversed upward from the bottom to the top wall 
for the S/D row spacing of 2.5 and 5. In both cases, 
the velocities measured in the central region of the 
enclosure showed no tendency to decrease as the probe 
approached either wall, falling within the standard 
deviation as calculated for that centerline velocity. 

Error in Directly Measured Variables 

The directly measured quantities are the axial 
location of the probe from the jet source and the 
nozzle diameter of the jet. The downstream position of 
the probe, X, is measured on a scale with a resolution 
of .0625 inches, while the diameter of the jet nozzle 
is presumed to be within .0001 inches. 

The Z and Y locations of the probe are never 
explicitly taken into account. It is assumed that the 
relative positioning of the probe due to the stepper 
operation of the traversing mechanism is accurate to 
within 5% of the last step, or 5% of 1.8 degrees. 1/200 
of a revolution of a 16 thread per inch shaft results 
in the linear travel of .0003 inches, with 5% of that 
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value being negligible. 

Error in Calculated Variables 

Fecause the hot wire probe is calibrated at a 
given temperature, it is not strictly a directly mea- 
sured variable. A 1 degree Fahrenheit change in the 
temperature from the 72 degree calibration temperature 
of the probe results in a .5% error in the anemometer 
reading. Given a 1 degree Fahrenheit resolution in 
reading the scale on an alcohol thermometer, the rela- 
tive error in the velocity reading is .5%. Coupled with 
.5% relative error associated with the bridge voltage 
calibration using a micromonometer to read wind tunnel 
flow speeds and a 2% relative error associated with the 
linearization of the bridge output, the total relative 
error for velocity readings is 3%. 

The internal error for Ui, or its standard devia- 
tion of the mean 

S.D. = [£{(Ui - U) 2 / 199 (200)}] 

is found to be 2%. At the jet source, the standard 
deviation is found to be less than 1.3%. 

The error in Q/Qo due to reading error is additive 
since the anemometer was used to determine each quanti- 
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ty. The external error alone then .if; 6%. The internal 
error in Q/Qo is also additive, for the same reason. 
For S/D = 2.5, 5 velocity distributions were determined 
at each section, giving an interna] error of 21/5 = 
.4%. This figure is added to 1.3% internal error found 
in the calculation of Co to give a total internal error 
of 1.7%. Added in quadrature to the external error, the 
total error in Q/Co for S/D = 2.5 is 6.3%. For S/D = 5, 
10, and 20 a minimum of 8 velocity distributions were 
established at each section, resulting in a total in- 
ternal error of 1.6% for Q/Qo. Again the internal and 
external errors are added in quadrature to yield a 
total error of 6.2% 

The entrainment formulas are dependent on X/D as 
well as Q/Cc, and at the near downstream sections of 
]6.4 and 32.8 the 3% and 1.5% relative errors in X/D 
must be added in quadrature to give the respective 
total errors of 7% and 6.3%, all other sections retain- 
ing the same value for the total error because of the 
diminishing relative error in X/D in dovnstrean sec- 


tions . 
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Velocity magnitude distributions 

The graphs presented in this appendix represent 
the velocity magnitude distributions from which the 
volume fluxes are calculated. Centerline velocity dis- 
tributions (Y/S = .5) are given for all sections, while 
velocity distributions taken at intermediate locations 
between adjacent jets, are incl uded for sections in 
which flow is still essentially three dimensional. 

The graphs plot, a dimensionless velocity (U/Uo) 
against a dimensionless distance from the jet center- 
line along the Z axis (Z/Z(l/2Um), where Z(l/2Um) de- 
notes the location at which, the velocity in a given 
section is half that of the centerline velocity. Urn, 
for the same section. 

The initial velocity profiles for each of the row 
of jets tested are else presented in Appendix B. The 
velocity is plotted against the no le radius, and as 
in the case of the total velocity, the distribution was 
recorded by traversing the hot wire probe (oriented 
along the Y axis) along the Z axis through the jet 
center! .ine . 
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FIGURE 8B1 Total velocity distribution for S/D 
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FIGURE 33Bi Total velocity distribution for S/D 
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FIGURE 47Bi Total velocity distribution for S/D 
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